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Abstract 

We consider the Maxwell-Lorentz equations, i.e., the equation of motion of a charged dust coupled to Maxwell's 
equations, on an arbitrary general-relativistic spacetime. We decompose this system of equations into evolution 
equations and constraints, and we demonstrate that the evolution equations are strongly hyperbolic. This 
result guarantees that the initial-value problem of the Maxwell-Lorentz equations is well-posed. We illustrate 
this general result with a discussion of spherically symmetric solutions on Minkowski spacetime. 

o ' 

CN ; 1 Introduction 

^^ • In this article we consider the equations of motion for a charged dust coupled to an electromagnetic field on a 

\my ', general-relativistic spacetime. These equations of motion are often refered to as the Maxwell-Lorentz equations, 

fT^ ' because they comprise the Maxwell equations for the electromagnetic field and the Lorentz force equation for the 

dust particles. Here, as usual, we use the word "dust" as synonymous to "pressure-less fluid" . As an alternative, 

-J ■ the name "cold fluid" is sometimes used. 

^h| The Maxwell-Lorentz equations have several interesting applications. E.g., they can be used for describing the 

electron component in a plasma. Among other things, this is of relevance for astrophysics, e.g., for modeling the 

bJQl Solar corona. In spite of the high temperature of the Solar corona, it is quite reasonable to model the electron 

component as a pressure-less (or "cold") fluid, because the density is low. The Maxwell-Lorentz equations have 

also been suggested for modeling a beam of charged particles in an accelerator, cf. Burton, Gratus and Tucker 

^ ' [I]- Special interest has been given to the Maxwell-Lorentz equations for rigid charge distributions on Minkowski 

CO I spacetime. The idea is to use this as a model for an extended charged classical particle, with the self-interaction 

•^ ■ taken into account. Interesting results have been found, e.g., by Spohn and collaborators [31 |3J |3] as well as by 

Bauer and Diirr [S]. An important goal of this line of research is to get some new insight into the point-particle 

limit and its notorious pathologies. 

r_^ I Several fundamental features of the Maxwell-Lorentz equations on Minkowski spacetime are discussed in Chapter 

5 of Parrott's text-book on electrodynamics [Bj. Among other things, Parrott discusses the initial- value problem for 

these equations. Using the Cauchy-Kowalevski theorem, he is able to demonstrate local existence and uniqueness 

of (analytic) solutions from analytic initial data. However, the assumption of analyticity is physically not quite 

satisfactory: Knowing an analytic function on an arbitrarily small neighborhood fixes the function on the whole 

j^ ' domain of convergence of the Taylor series, whereas physical measurements in an arbitrarily small neighborhood 

H I cannot determine any results outside of the causal future of this neighborhood. Also, it is often convenient to 

consider fields which are zero on an open subset U and non-zero on the complement of U; such fields cannot be 

analytic near the boundary oiU. For these reasons it is desirable to have an existence and uniqueness theorem for 

non-analytic data, as noted by Parrott himself. It is the main purpose of this paper to establish such a result. 

To that end we decompose the Maxwell-Lorentz equations, on an arbitrary general-relativistic spacetime, into 
evolution equations and constraints. The evolution equations can be written as a quasi-linear system of first order 
partial differential equations. We derive the characteristic equation (dispersion relation), and we demonstrate that 
the system of evolution equations is strongly hyperbolic, which is a necessary and sufficient condition for the initial- 
value problem to be well-posed. This means that a solution locally exists, that it is unique, and that it depends 
continuously on the data in an appropriate topology. This result is of relevance for solving the Maxwell-Lorentz 
equations numerically: If a solution does not exist it is pointless to write a numerical code that tries to approximate 
it; if the solution is not unique the code does not know which solution to approximate; and if the dependence on 
the data is discontinuous it is impossible to estimate the error of a numerical solution. 

After establishing our general existence and uniqueness result, we illustrate it with a discussion of spherically 
symmetric solutions on Minkowski spacetime. 
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2 The Maxwell-Lorentz equations 

We consider the coupled equations of motion of a charged dust and an electromagnetic field on a general-relativistic 
spacetime, i.e., on a 4-dimensional manifold M with a Lorentzian metric g. We allow for an external electric current 
density J^, and we assume that the dust consists of only one species of particles, with charge q and mass m. In the 
following we assume that the second-rank tensor field g = gijdx^ (E}dx^ , the one-form J'^ = Jfdx^ and the constants 
q and m are given. We use standard index notation with Einstein's summation convention for latin indices taking 
values 0, 1, 2, 3 and for greek indices taking values 1, 2, 3. As usual, the contravariant metric components are defined 
by g'^^gij = 6j. The metric has signature (— , -I-, -I-, -I-) and its Levi-Civita connection is denoted V. We use units 
making Eq and fiQ (and, thus, the vacuum speed of light c) equal to 1. 

The equations of motion consist of Maxwell's equations for the electromagnetic field and the Lorentz force 
equation for the dust particles, 

d^Fjk + ^JFk^ + dkF,, = , (1) 

g''y^F,k = Jk - qngkiU' , (2) 

W\/,W ^ ^g^'^FkiU' . (3) 

m 

We refer to these equations as to the Maxwell-Lorentz equations. The field variables (i.e., the "unknowns" in these 
equations) are 

• the electromagnetic field strength Fijdx'^®dx^ which is an antisymmetric second-rank tensor field, Fij — —Fji ; 

• the four-velocity V^di of the dust which is a vector field, normalised to gijWU^ ~ —1 . 

• the (proper) number density n of the dust which is a non-negative scalar function; 

The electromagnetic field is to be interpreted as the field produced by the external current plus the field produced 
by the charged dust plus, possibly, an additional source-free field. 

If n = on an open subset of the spacetime, C/* is physically undetermined on this subset. In this case, the 
electromagnetic field is determined by ((ij and ([2]), i.e., by Maxwell's equations with a given source J^. As the 
initial- value problem for this system of equations is well-known, we assume in the following that n > on our 
spacetime manifold. (Or, to put this another way, we restrict our spacetime to an open subset on which n > 0.) 

Note that ([3]) implies U^di[gjkU^U^) = 0. This guarantees that, if the condition gjkUW' = —1 is satisfied 
at one point, then it is satisfied along the whole integral curve of U^di that passes through this point. A similar 
statement is true for the condition n > 0, provided that the external current is conserved, g^^^iJ^ = 0- To prove 
this, we read from ([2]) that the conservation law for the external current implies Vfc(n[/ ) = and thus 

U^dun + nykU^ = Q . (4) 

If x{t) is an integral curve of the vector field U^di, integration of ^ yields 

n{x{t2)) = n{x{ti)) exp{ - / ' (VfeC/^') {x{t)) dt\ (5) 

so n cannot change sign along an integral curve of U^di. This is an important observation in view of the initial value 
problem. It makes sure that, for a solution of the Maxwell-Lorentz equations, it suffices to check the conditions 
gijWU^ — —1 and n > on the initial hypersurface. 

To illustrate the physical content of the Maxwell-Lorentz equations we briefly mention two different applications. 
Firstly, we can use these equations for modeling plasma waves in the Solar corona. In this case we would choose for 
(M, g) the Schwarzschild spacetime, for J'^ the electric current density of the ion component which could be assumed 
static, and for q and m the charge and the mass of the electron. The fleld variables are the total electromagnetic 
field, the four- velocity of the electrons and the number density of the electrons. Secondly, we can consider an 
electron beam in an accelerator. Then we would choose for (M, g) an open subset of the Minkowski spacetime (to 
model the spacetime region inside the accelerator pipe), for J*^ the zero one- form (as there is no external current 
inside the accelerator pipe), and for q and m the charge and the mass of the electron. The field variables are again 
the total electromagnetic field (i.e., the field produced by the dust and the accelerating field which has no sources 
inside the pipe), the four- velocity of the electrons and the number density of the electrons. 



3 Decomposition of the Maxwell-Lorentz equations into evolution equa- 
tions and constraints 

We fix a spacetinie point P and we choose local coordinates, on an appropriate neighborhood of P, such that the 
hypersurfaces x'^ = constant are spacelike and the x*^— lines are perpendicular to these hypersurfaces. We want to 
discuss if the prescription of appropriate initial data on the hypersurface a;° — constant through P determines a 
unique solution of the Maxwell-Lorentz equations on some neighborhood of P. To that end we have to decom- 
pose the Maxwell-Lorentz equations into evolution equations (i.e., equations that do contain 9o derivatives) and 
constraints (i.e., equations that do not contain 9o derivatives). 

There are two constraints, given by equation ([l} with all three indices spatial and by equation ^ with k = Q. 
After expressing covariant derivatives in terms of partial derivatives, introducing the ChristofFel symbols of the 
metric tensor, these two equations take the following form. 

a^i^.a + d,F„^ + d,F^, = , (6) 

9^"- {d^F^o - r;,F^o - r^oF,, - r%F,o) ^ 4-qn gooU'' ■ (7) 

The remaining equations are the evolution equations: Evaluating equation ([T|) with one of the indices equal to 
gives three scalar equations, 

doF^, + d^F^o + d^Fo^ = . (8) 

Evaluating equation ^ with a spatial index, k — a, gives another three scalar equations. After eliminating the 
number density n, with the help of ([7]), these equations read 

JJP je JJP 

~lp9<yp9'^''{^J,i.FTO + ^%F^P + ^fj,oP'^o) + roo-Fo<T + r^JpF^^ + FJ^^Fom + '^ ~ Jp dap Jq ■ (9) 



Finally, evaluating equation ([3]) with a spatial index, j — a, gives us again three scalar equations. With the help 
of equation ([3]) with j = 0, these equations can be rewritten as 

jja JJP jja- jjn JJf'JJ'' 

Pi I ft _ _ po- _ 9pcr ^ _ per ^ ^ 

U° U° "U^ ^ "° ^^opjjo ^ P'^ /ijo\2 

Equation ^ with j = can be ignored because it is automatically satisfied if the other equations are, owing to 
the normalisation condition gijWU^ — —1 . So the system of evolution equations is given by nine scalar equations, 
(ISI), (0) and ([TU)) . Note that the number density n is not a dynamical variable; n has been eliminated from the 
evolution equations and is algebraically determined by the other field variables via (jT]). 

To write the evolution equations in a more convenient form, we decompose the electromagnetic field into its 
electric and magnetic components, E^ and -B", with respect to the chosen coordinate system, 

F,o = E^ and Frp = ErpuB" . (11) 

Here Sa-pu is the three-dimensional covariant e-symbol, defined by being totally antisymmetric and satisfying 
^123 = 1. Moreover, we introduce the 3-velocity v"' of the dust, 

Wd, ^ C/O {do + v'^d,) , (12) 

Note that then C/'^ is determined by the normalisation condition (f7")^(.goo + v^^v'^ g^^) = — 1 . If we use the three- 
dimensional contravariant e-symbol e'^'"', defined by being totally antisymmetric and satisfying e^'^^ — 1, we can 



write the evolution equations in the following final form. (We use the identity e^'^'^e^^r = 2(5^-) 

doB" + e^'^-'df.E, = , (13) 

+ vP g^pg^-^ {Tl^Er + T^o^.pAi?^ + V%E,) - ^ + v' 9^p Jo + ^ooEa - Ko^,..B'' + T^^E^ , (14) 

+ ^ v/-5oo - vt^v-g^,, [g^^ ( Er + Erpxy^B^ ) + g°%''v^Er } . (15) 

With the same notation, the constraints ([6]) and (O read 

dpBf^ = , (16) 

- .900 ^ ^ 

Note that (|13p implies dodaB'^ = 0. Hence a solution of the evolution equations satisfies P^ everywhere if it 

satisfies P^ on the initial hypersurface x" = constant. As a consequence, the solutions to the Maxwell-Lorentz 

equations can be found by solving the initial- value problem of the evolution equations, with initial values for Ea-, 

B^ and v^ that satisfy (flB)) and that give a non- negative number density n via ((T7)) . Recall that, by ([5]), it suffices 

to check the condition of n being positive on the initial hypersurface. 

If we write 

^E \ IB^^ 

E=\E^\ , B=\B^\ , v^\v'\ ^ (18) 

?3. 



.E.} \B^ 



the evolution equations (J13p , (J14p and (jlSp can be written in matrix form as a differential equation for a 9-component 
column vector, 

do Ib] - LPOpIb] = I ml] . (19) 

Here, L , L , L are 9x9 matrices and mi, TO2, rn^ are 3-componcnt column vectors. All of them depend on the 
spacetime coordinates, via the given background fields. In addition, the L'^ depend on v and the rhp depend on 
all dynamical variables {E,B,v), but not on their derivatives. Hence, P^ is a quasi-linear system of first-order 
differential equations for {E,B,v). 

We can choose our coordinates such that, at the fixed spacetime point P, we have (g^jy) = diag(— 1, 1, 1, 1) and 
rj^^ = 0. Then we can read from dni), (Ull) and ^5^ that, at this point. 



LP = \ AP (20) 





where 




e- = 1 , e^ ^ , (21) 



(22) 





The fhp are given, at the point P, by 



mi == /"= + J^v , m2 = , ma = — ^1 - v^v>^6rx [E + s.^^v^'B^e" ~ v^Epv) . (23) 

If (M, 5) is Minkowski spacetime, and if we use standard Minkowski coordinates, ((20)) and ((23)) hold not only at 
one point but everywhere. 

4 Strong hyperbolicity of the Maxwell-Lorentz equations 

In this section we will demonstrate that the initial- value problem for the evolution equations of the Maxwell-Lorentz 
system is well-posed. If we assume that the background fields are analytic, and that we give analytic initial data 
for (E,B,v), local existence and uniqueness of a solution follows from the Cauchy-Kowalevsky theorem; this was 
the result found by Parrott [6 . (Parrott only considered the case that {M,g) is Minkowski spacetime and that 
J"^ vanishes, but the result carries over immediately to our more general case.) For the reasons outlined in the 
introduction, we want to have an existence and uniqueness statement for non-analytic data. To that end we have 
to recall some terminology and some results from the theory of partial differential equations. 

For a quasi-linear system of first-order partial differential equations of the form ((TS)) . the left-hand side of (H^j) 
is called the principal part. Existence and uniqueness of solutions is completely determined by the principal part 
and hence by the matrices L^. ((TS)) is called hyperbolic if for any p = PpC^ G M.^ the eigenvalues of the 9x9 
matrix PpL'' are real. (|19p is called strongly hyperbolic or symmetrisable if, in addition, for each eigenvalue the 
geometric multiplicity equals the algebraic multiplicity, i.e., iippL^ has nine linearly independent real eigenvectors. 
An equivalent condition is that there exists a matrix H{p) such that H{p)ppL'' is symmetric. 

According to a fundamental theorem from the theory of partial differential equations (see, e.g. Taylor [7], 
Theorem 5.2.D), strong hyperbolicity guarantees that the initial- value problem is well-posed for Sobolev H^ data. 
(A function is of class iJ" if its s**^ derivative exists almost everywhere and is locally square-integrable.) In general, 
s must be chosen greater than {n + 2)/2, where n is the dimension of the initial hypersurface. If adapted to our 
case, the theorem can be stated in the following way. 

Assume that (|19p is strongly hyperbolic, with the matrices L'' and the vectors rfip smoothly depending on the 
spacetime coordinates and on the dynamical variables {E,B,v). Choose initial data of class iP* for {E,B,v) on 
the hypersurface x^ ^constant through the spacetime point P, where s > 3. Then 

• a solution to ((TS)) with the prescribed initial data exists on a neighborhood of P ; 

• on this neighborhood, the solution is uniquely determined by the initial data; 

• the solution depends continuously, with respect to the iP* topology, on the initial data. 

It is now our goal to demonstrate that the evolution equations of the Maxwell-Lorentz system are strongly 
hyperbolic. If this result has been established, we know that the initial-value problem is well-posed, provided that 
our background fields g^^ and Jfi depend smoothly on the spacetime coordinates and the initial data are of class 
H" with s > 3. 

As strong hyperbolicity is an algebraic property that is to be checked pointwise, we may assume that the 
matrices L'' are given by ((^H)) . hence 




PpLP = PpAP ■ (24) 



We have to show that this matrix has nine linearly independent real eigenvectors for each p = Ppe^ S M'^. We first 
determine the eigenvalues, po is an eigenvalue of PpL'' if the characteristic equation 

Aei{pol~PpLP) = (25) 

is satisfied. If this equation holds, pidx^ = p^dx^ + Pp,dx^ is called a characteristic covector of the evolution 
equations. Upon inserting PpL^ from ((24p into ((25p . the determinant can be calculated with the help of some 
elementary algebra. As a result of this calculation, the characteristic equation reads 

{pI ~ S^'Pt^Puf (Po + PpVf po = , (26) 



so the eigenvalues are y/Sf^'^p^p^, —^/^^''PiJ.Pv, —PpV^ and 0. As all eigenvalues are real, this demonstrates that 
the evolution equations are hyperbolic. 

We will now prove that they are even strongly hyperbolic. To that end we have to determine the corresponding 
eigenvectors. As the case p = is trivial, we may assume that Py^O. Owing to the homogeneity of the eigenvalue 
equation it suffices to consider the case that S^'^p^j^p^ = 1. It is then straight-forward to verify that 




■ = UpcP e W^ , S^'^Uf.p^ = 



(27) 



is a two-dimensional eigenspace of the eigenvalue pq — 1, 



£2 = 




= wpcp e m.^ , S'^'wap^ = 



(28) 



is a two-dimensional eigenspace of the eigenvalue pa = —1, 

PpvPy 
^3 == "^ I -apnAPy 



-appj 

z 



a G M, z e 



(29) 



is a four-dimensional eigenspace of the eigenvalue po — ~Ppv'' , and 



/3e 




(30) 



The eigenspaces f 1 , E2, f 3 and £4 span the whole space, 
= or ppv'' = ±1). This demonstrates that the evolution 



is a one-dimensional eigenspace of the eigenvalue po = 
even in the case of further degeneracy (i.e., even iippV' 
equations are strongly hyperbolic. 

This concludes the local existence-and-uniqueness proof for solutions {E,B,v) to the evolution equations. For 
solving the full set of Maxwell-Lorentz equations, we have to single out those solutions to the evolution equations 
that satisfy the constraints. As a consequence, not all eigenvectors in £1, £2, £3 and £4 are associated with solutions 
of the full set of Maxwell-Lorentz equations. The condition n > does not restrict the eigenvectors. However, 
equation p^ reduces the eigenspaces according to 




^£a 



P^B^ = 



A = l,2,3,4. 



(31) 



By inspection we find that £1 — £1, £2 — £2, £3 = £3 whereas £4 consists of the zero- vector only, so po = is 
no longer an eigenvalue if the constraints are taken into account. This reduces the characteristic equation of the 
evolution equations (j26p to the characteristic equation of the full set of Maxwell-Lorentz equations. 



(Po - S'^'^P^P.f iPo 



v'PpY 







In physical terms, (j32p gives us the three branches of the dispersion relation, 

Po = ^/S^'^PpPy , po ^ -\/5'^''PpP,,, PQ = -vPpp 



(32) 



(33) 



From the corresponding eigenspaces we read that the first two branches are associated with the (past-oriented and 
future-oriented, respectively) propagation of the electromagnetic field, whereas the third branch is associated with 
the propagation of the dust. If we calculate the group velocity (or ray velocity) 



9(-Po) 
dp,. 



(34) 



for each branch of the dispersion relation we see that the electromagnetic field can propagate in any spatial direction 
with the vacuum speed of light ( v^'^Vgj.^^^ = 1), whereas the dust propagates with v^j. = w^. 



5 Example: Spherically symmetric solutions on Minkowski spacetime 

Knowing that the initial-value problem admits a miique solution is one thing, actually determining this solution 
is another. In most cases one can find the solution only numerically. However, in some cases of high symmetry 
it is actually possible to work out some features of the solution explicitly. Parrot [SI considers two such cases: (i) 
spherically symmetric solutions on Minkowski spacetime, and (ii) solutions on Minkowski spacetime that depend 
only on one of the three spatial Cartesian coordinates {x,y,z). The first case is concerned with expanding or 
collapsing shells of charge, and the second case with moving "walls" of charge. Burton, Gratus and Tucker also 
considered the first case [8, and the second case 1^. 

Here we want to revisit spherically symmetric solutions. Using Parrott's treatment of this case as a general 
guide, we will get some results that go beyond his findings. In particular, we will solve the differential equations 
for the flow lines of the dust explicitly, and we will get some additional insight into the qualitative features of the 
solutions. Although spherically symmetric fields are somewhat over-idealised, we find it instructive to study this 
case. As we will see, it illustrates the way in which the constraints have to be built in and it demonstrates how 
caustics are formed, thereby preventing the solution from being global. 

We assume that the background metric is the Minkowski metric and that the external current vanishes, J^ = 0. 
In spherical polar coordinates (x'^ , x^ , x'^ , x^) = (t,r,'d,if), the metric reads 

gijdx'dx^ = -df + dr^ + r^ {siT?ddip^ + dd^) . (35) 

For spherically symmetric situations, the constraints equation p6p can hold with a non-zero magnetic field only if 
there is a magnetic monopole at r = 0. We exclude the existence of magnetic monopoles and, thus, assume that 
B'^ = 0. We are then left with only two dynamical variables, 

Ef,dx^' = E{r,t)dr , (36) 

v^d^ = v{r,t)dr, (37) 

and the evolution equations p3)). (fT4)) and (fT5| are identically satisfied, apart from the following two scalar 
equations, 

dE(r.t) , , dEir.t) 2 , ^ ^, , 

\, ' +v{r,t) 1' ' = --v{r,t)E{r,t), 38 

ot or r 

dv{r,t) dv{r,t) q 



n, + <^^ t) ^^^ = - Vi-vir,t)^ Eir, t) , (39) 

ot or m 



The constraints (TTSl) and dTT]) reduce to a single scalar equation, 

q n{r, t) = ^l - v{r, t^ ( ^^^ + ^ E{r, t) ) , (40) 

which gives the charge density q n{r, t) in terms of the dynamical variables E{r, t) and u(r, t). Note that the electric 
field is to be interpreted as the field produced by the charged dust plus, possibly, a source-free field. The only 
spherically symmetric electric field that is source-free, on a spherical shell ri < r < r2, is a Coulomb field; i.e., it 
can be realised by a point charge at r = 0. 

To solve the system of equations ([55]) . ([M]). and PO]) we would have to give initial values E{r,to) and v{r,to) 
on a certain interval ri < r < r2 . With such initial data given, our existence and uniqueness result guarantees that 
there is a unique local solution to ([55)) and ([5^ . From all possible initial values we would then have to discard all 
those for which (|40p gives a negative number density n(r, t) on the initial hypersurface. 

The latter constraint is a bit awkward. It would be more convenient if we could give the initial data in terms of 
n{r,to) directly. Also, if we want to explicitly integrate ((38)) and (|39)) we have to face the problem that these two 
equations are coupled. It would be nicer to have an equation for either the electric field or the flowlines of the dust 
alone. Both goals can, indeed, be achieved if we rearrange our system of equations appropriately. To that end we 
first define 

Q{r,t) = Anr^ E{r,t) . (41) 

This definition follows Parrott, adapted to our choice of units. (Parrott has a factor An in Maxwell's equations, we 
don't.) Then (j^O)) can be rewritten as 



,„,,„.^IIlMiMid)^ (42) 

4 TT r^ or 



Upon integration, this yields 

n{s, t) s^ ds 



Q{r,t) = Qiro,t) + 4^q / ^ ' (43) 

ro V -*- ~ ^(■Sj^j 



where tq is any chosen radius value. As qn{r,t) is the proper charge density, q n{r, t)/ y^l — u(r, t)^ is the charge 
density in the chosen inertial frame. Hence we read from equation (j43p that Q{r, t) — Q{rQ, t) is the total charge in 
the shell between radius rg and radius r at time t. After replacing the original pair of dynamical variables [E, v) 
by (Q,w), the evolution equations ((38|) and (p9)) read 

2^+„(..„5O&i).0, (44, 



^ + „,,, ,) ^ = ± ^' -;^'f Q(r. t). (45) 

at or m 47rr^ 

If we evaluate these equations along a flow line [r{t),(p{t) = •^Q,-&{t) — -do) of the dust, such that v[r{t),t) = 
dr{t)/dt, we find 

dQ(r(t),t) , , 

^l^ ' = , (46) 

dt m 4 TT r^ 

The first equation says that Q is conserved along each flow line of the dust. With this information at hand, the 
second equation can be viewed as a second-order differential equation for the flow lines alone. We can, thus, solve 
our system of equations in the following way. 

Fix a time io Sind a radius value tq. Prescribe, on an interval ri < r < r2, initial values 

n(r,io) = nn{r) , no(r) > , (48) 

v{r,to) = voir) , vo{rf < 1, (49) 

E{ro,to) = ^0- (50) 

The freedom of choosing a value for Eq corresponds to the freedom of adding a Coulomb field, i.e., a point charge 
at r = 0. 

The data (|^5|) . (P9|) and (|50p determine initial values for the function Q, namely 

Q(r,io) = ^TrrlEo + 4nq f ""f'^'^f' . (51) 

Jr„ Vl-'yo(s)^ 

Then (|T7)) gives us a second-order differential equation for the flow lines of the dust, 

d^rjt) _ ^ r , /dr(t)\2 ^3/2 Q(/,to) 
dt2 ^ m I Wt J i r(i)2 ' ^ ^ 

that has to be solved with initial conditions 

dv 
r(io) = r' and —{to) = vo{r') (53) 

for all ri < r' < r2. If we have determined the flow lines, our dynamical problem is essentially solved. It is true 
that knowing the flow lines does not allow us to write closed form expressions for the other dynamical variables. 
However, all the relevant information is then known: The function Q, which determines the charge in any spherical 
shell, is fixed by being constant along the fiow lines; thereupon the electric field is fixed by (HT|) . The solution exists 
and is unique until the flow lines begin crossing over and forming caustics. 

So our problem has been reduced, in essence, to solving the differential equation (j52p . This can actually be 
done explicitly. After multiplying both sides with dr{t)/dt, ([5^ can be rewritten as 

d/r, {dr{t)Y\-'/^\ qQjr'M) d( 1 x 



This gives us a first integral, 



I \ dt J } mr{t) ^ 

The value of the constant of motion A(r') is determined by evaluating (|55p at t — to, 



mr 



(55) 



(56) 



Equation (j55p has the interesting consequence that a dust particle can reach the origin, r{t) = 0, only if it arrives 
there at the speed of light, i^dr{t)/dt) = 1. The only exception are trajectories with Q{r\to) = which, by 
([5^ . are straight lines. Note that trajectories with qQ{r',to) > are accelerating outwards and trajectories with 
qQ{r\to) < are accelerating inwards. 

The first-order differential equation ([55]) can be solved by separation of variables. We find 



qQ{r',h) 



±*= ^^^^^^^[^VW)\IZir')r^ 



2 ^Q^''^''K \{r')r I g'Q("^^°)' 



m 



(57) 



■2Z{r')r 



^qQJr'M^^^, 



B{r') 



\iZ{r') = A{r'f-l>Q, 



± t = 



qQjr'M) 
m ( - Z{r')) 
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if Z{r') = A{r'f - 1 < 0, and 



±t^M-2^i^^^^^Air')r^l^^'^ ' 



iA{r') 3qQir\to) 



B{r') 



(59) 



if Z{r') = A{r'Y — 1=0. In all three cases, the integration constant B{r') is determined by setting t ^ to and 
r = r' , and the sign on the left-hand side must be chosen positive for outward going and negative for inward going 
parts of the trajectory. 

For any specific choice of the initial conditions, ([57]). ([55)1 and (l5^ give us the fiow lines of the dust which 
can then be plotted in an (r, t) spacetime diagram. For the sake of illustration, we choose the following initial 
conditions on the interval < r < ro, 



n{r,to) = a{ro - rf 

vir,to) = , 

E{ro,to) = Eo, 



(60) 
(61) 
(62) 



where a is a positive constant and Eq is a constant that may be chosen positive, negative, or zero. In this case 
our system describes a charged ball of radius ro that is initially at rest; in addition, there is a point charge at the 
origin whose value depends on the chosen value of Eo- 

Figures [1] [5] and [3] show Maple plots of the flow lines in an (r, t) spacetime diagram for three qualitatively 
different cases, depending on the chosen value of Eo- As discussed in the captions of the figures, this example 
nicely demonstrates the role of the point charge at the origin (and, thereby, of the fact that in general the choice 
of initial values for the Maxwell-Lorentz equations includes the choice of a source- free electromagnetic field). It 
also illustrates how a caustic is formed (which is the main reason why, in general, existence and uniqueness of a 
solution is guaranteed only locally near the initial hypersurface) . 




Figure 1: If qEq is negative, qQ(r,to) is negative on the whole interval < r < tq. In this case there is a point 
charge at the centre whose charge has the opposite sign as the charge of the dust. The Coulomb attraction of this 
point charge is so strong that it overcomes the Coulomb repulsion of the dust particles and sucks the whole ball of 
charge into the origin in a finite time. All dust particles arrive at the origin with the speed of light. Although the 
initial condition for the dust's density, n{r,to), admits a regular extension into the origin, the density is singular at 
the origin for all later times. If one assumes that the dust particles can pass through each other when they meet 
at the origin, the ball would re-expand after its collapse. 
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Figure 2: If < qEq < g^arg/SO, there is a radius value r between and tq for which Q{r,to) becomes zero. At 
this radius, the ball is stationary. The outer part expands because here the Coulomb repulsion of the dust particles 
is stronger than the Coulomb attraction of the point charge towards the origin; the inner part collapses because 
here it is vice versa. Again, any dust particle that reaches the origin arrives there at the speed of light, and the 
density n(r, t) is singular at r = for all t > to- 



ll 




Figure 3: If qEo > (/^arg/SO, the point eharge at the centre has the same sign as the charge of the dust, so it 
enhances the tendency of the dust to expand outwards. The inner parts of the ball move faster than the outer 
parts and overtake them, thereby forming a caustic. The solution to the Maxwell-Lorentz equations breaks down 
when flow lines cross over. As guaranteed by our general theorem, the solution is well-defined and unique on a 
neighborhood of the initial hypersurface. In this case, the latter is a snapshot of a punctured ball. 



12 



6 Conclusion 

In this article we have demonstrated that the initial- value problem of the Maxwell-Lorentz equations is well-posed, 
and we have illustrated some qualitative features of these equations by studying spherically symmetric solutions on 
Minkowski spacetime. In view of applications, it is desirable to generalise our result in two directions. Firstly, one 
would like to include a pressure term into the equations, thereby allowing for situations where temperature and 
density are so high that the model of a "dust" (or "cold fluid") is no longer applicable. Secondly, one would like to 
consider mixtures of two or more fluids which are dynamically coupled together. This would allow for considering 
two-fluid models of plasmas, treating both the electron component and the ion component as dynamical fields. In 
particular the first generalisation would require a considerable change in the mathematical techniques. Whereas 
for our treatment of a dust it was crucial that the density could be eliminated from the evolution equations, this 
would no longer be true for a fiuid with pressure. One would have to use either the density or the pressure as an 
additional dynamical variable, and to eliminate the other one with the help of an equation of state. One would also 
expect that a pressure term would change the behaviour of a collapsing ball of charge. In particular, the somewhat 
pathological behaviour near the centre, as discussed in Chapter [3 is likely to be modified by a pressure term. 

Acknowledgment 

VP wishes to thank Nico Giulini for clarifying comments on the interpretation of solutions to the Maxwell-Lorentz 
equations. VP has also profited from discussions on the subject with his Lancaster colleagues, in particular with 
David Burton, Jonathan Gratus and Robin Tucker. 

References 

[1] D. A. Burton, J. Gratus, and R. W. Tucker. Asymptotic analysis of ultra-relativistic charge. Ann. Phys. (NY), 
322:599-630, 2007. 

[2] H. Spohn. Dynamics of charged particles and their radiation field. Cambridge University Press, Cambridge, 
2004. 

[3] A. Komech and H. Spohn. Long-time asymptotics for the coupled Maxwell-Lorentz equations. Commun. Partial 
Diff. Equat., 25:559-584, 2000. 

[4] M. Kunze and H. Spohn. Adiabatic limit for the Maxwell-Lorentz equations. Ann. H. Poincare, 1:625-654, 
2000. 

[5] G. Bauer and D. Diirr. The Maxwell-Lorentz system of a rigid charge. Ann. H. Poincare, 2:179-196, 2001. 

[6] S. Parrott. Relativistic electrodynamics and differential geometry. Springer, New York, 1987. 

[7] M. E. Taylor. Pseudodifferential operators and nonlinear PDE. Birkhauser, Boston, 1991. 

[8] D. A. Burton, J. Gratus, and R. W. Tucker. Multiple currents in charged continua. J. Phys. A, 40:811-829, 
2007. 



13 



